# Computes squared losses 
weighted.squared.loss <-function(pin){ 
	ein <- A.mat[q.pos,z.pos]%*%pin - P[q.pos,t]
	#ein <- A.mat%*%pin - p
	#wgt <- diag(length(q.pos))
	return(as.vector(t(ein)%*%diag(W[q.pos,t])%*%ein))
	#return(as.vector(t(ein)%*%wgt%*%ein))
}